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Abstract 

The axial N — A(1232) transition form factors are calculated within the light cone 
QCD sum rules method. A comparison of our results with the predictions of lattice 
theory and quark model is presented. 
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1 Introduction 

Form factors have great importance in the investigation of the internal struc- 
ture of baryons. The inner structure of the nucleon is encoded in several form 
factors. The electromagnetic form factors of the nucleon, which parameterize 
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the matrix element of the electromagnetic current operators, are measured in 
a wide range of values [T1I21I3] . The interest in nucleon form factor has been 
renewed by recent progress in experimental physics, where it became possible 
to get polarized beams and polarized targets. This possibility opened up two 
new methods [1] to measure the ratio of the electric and magnetic form factors, 
namely polarization transfer and beam-target asymmetry. New experimental 
data obtained from e + p ^ e + p reaction, which is performed at JLAB [5], 
gave the result that the ratio ij^pGeIGm decreases from unity substantially for 
high values. Not only form factors relevant for the diagonal transitions, 
but also the electromagnetic transition form factors for electro-production of 
the A has also been the subject of recent experimental [5|6|7] and theoretical 
studies [B]. 

The main advantage of the nucleon to A transition is that, the A is a dominant 
nucleon resonance and its identification is easy since the spin parity selection 
rules provide us information about the wave function of these baryons. The 
N — A transition form factors due to the weak axial current can give valu- 
able information about the structure of the baryons, complementary to that 
obtained from electromagnetic transition. For example, measurement of axial 
form factors for N — A transition, allows us to check off diagonal Goldberger- 
Treiman relation, order of conservation of the axial current, etc. Weak axial 
N — A transition form factors are investigated in neutrino (or charged lep- 
ton) scattering on deuterium or hydrogen in the A region experiments. New 
information on the weak axial form factors is expected from parity-violating 
electron scattering experiments planned at Jefferson Laboratory [5]. 

In the present work, we calculate the axial N A transition form factors 
in light cone QCD sum rules (LCQSR). In the light cone QCD sum rules 
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method, operator product expansion is carried out near the hght cone ~ 0, 
and non perturbative dynamics is parameterized by the hght cone distribu- 
tion amphtudes that determine the matrix elements of the nonlocal operators 
between vacuum and one particle states. The expansion near the light cone 
is an expansion in the twist of the operators rather than the dimension as in 
the traditional QCD sum rules (for more about this method see e.g. PHIITT] 
and references therein). The light cone distribution amplitudes of the proton 
is calculated in pHl25f26] . Using the light cone distribution amplitudes of the 
proton semileptonic Af, — > pii/ [15] decay, scalar form factor of the proton 
[T6IT7] . axial and induced pseudo scalar form factors of the nucleon [TSfTQ] . 
H — n form factors [20j are studied. In [21j, the distribution amplitudes of A 
in leading conformal spin is calculated and the obtained amplitudes are used 
to study the Aiu decay. Note that these form factors are calculated in 

lattice QCD in p^][T3] . and in chiral constituent quark model in [Ti] . 

The plan of this work is as follows: In section 2 we consider the generic correla- 
tor function and present the LCQSR formalism. In this section we obtain the 
LCQSR for the axial to A transition form factors. The numerical analysis 
and discussion is presented in section 3. 



2 Sum Rules for the Axial Nucleon to Delta Transition Form Fac- 
tors 

In the present section, light cone QCD sum rules for the axial A^ to A transition 
form factors are derived. These form factors are defined by the matrix element 
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of the axial current (in our case isovector part of the axial current), i.e. 
- 

Jl[x) = ij{x)-f^-f5—ij{x) (1) 

between the nucleon state with momentum p and the A-state with momen- 
tum p' = p — q, where is the third Pauli matrix, and q is the transferred 
momentum. Hence, the axial — > A transition form factors are defined by 
the matrix element < A{p',s') \ \ N{p,s) >. This to A weak matrix 
element can be expressed in terms of four invariant transition form factors as 
follows [22i.23j : 



< A{p', S') I 4 I N{p, S)>= IU\P\ S'){{^^^1, + 9±^p'^)i^g^^g^^ 



C^iq^) 



m 



N 



(2) 



where u^{p\ s') is the Rarita-Schwinger spinor for A and u{p, s) is the nucleon 
spinor. Note that partial conservation of axial current (PCAC) and pion meson 
dominance leads to the relation 
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For the calculation of the above mentioned form factors within the light cone 
QCD sum rule method we consider the matrix element in which one of hadrons 
is described by an interpolating current with quantum numbers of this hadron 
and another one is represented by the state vector. For the construction of the 
light cone sum rules, the distribution amplitudes (DA's) of state vector hadron 
is needed. The nucleon distribution amplitudes for all three quark operators 
are calculated in p^25f26] and DA's for A isobar is not yet calculated. For 



4 



this reason, we consider the following correlator function where nucleon is 
described by state vector 

n;..(p, q) = tj rfW'?" < I T{r]^{0)4{x)} \ N{p) > (4) 

The interpolating current for the A"*" isobar is chosen in the form [27] 

rj.iO) = -^e'^"^[2iu'^^iO)C^,d\0)X{0) + K^(0)C7y (0))rf^(0)] (5) 

and the axial current is: 

Juix) = ^[M(a;)7^75M(x) -d{x)'y^'y5d{x)] (6) 

In order to construct the sum rules for the form factors, the correlation func- 
tion will be represented in two different forms, i.e. in terms of hadron pa- 
rameters and in terms of the quark-gluon parameters. Let us first consider 
the representation of the correlation function in terms of hadron parameters 
(phenomenological part). Phenomenological part of the correlation function 
can be obtained by inserting the complete set of hadrons with the same quan- 
tum numbers of ^7^(0) inside the correlation function. Saturating the correlator 
function with these hadrons and isolating the contribution coming from ground 
state hadron we get 

. j: <o\n,.\^nM><AYs')\J,\N(p.s)> ^ , 

where ttt-a is A"*" mass and ■ • • represents contributions from higher states and 
the continuum. The matrix element < A+ | J;^(x) | > entering in Eq. ([7]) is 
given by Eq. ([2]) and the matrix element < | 77^(0) | A+ > is determined as 

<0\r],{0)\A+>=X^u^{p',s') (8) 
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where Aa is the residue of A+ isobar and the value of the residue Aa is deter- 



mined from the two point sum rules [27|28II29|30|31] . Performing summation 
over spin of Rarita-Schwinger spinor using the formula 



u^^ip', s)u^{p', s) 



+ 



}(9) 



and using Eqs. (E]), ([7]), (IE]) and ([9]) the contribution of to the correlation 
function can be written as 



—iX/ 



771/ 



p 



/2 



(/ + TTIa) 



2p>A ^ p;,7a - Pa7/. 



3mA 



rriN 



-7a + 



m 



-Pc 



N 



+ 



(10) 



From Eq. ( iTOl) it follows that the correlator function contains numerous ten- 
sor structures and not all of them are independent. The dependence can be 
removed by ordering gamma matrices in a specific order. In this work the 
ordering 7^ ^''-f,, /j is chosen. With this ordering the correlation function 
becomes 



^f^uip,q) 



K,. 4ci{q') + g,. a + gM + ^a) 



m^ — 



m^v 



-q^ip + m/^)p^ 2 ^ 



m 



N 



m 



N 



+ other structures with 7^ at the beginning or which are 

proportional to p^. (11) 



The reason why the structures ~ p'^ and structures with 7^ at the beginning 
is not considered is follows. The interpolating current r^^ couples not only to 
spin-parity (3/2)+ states, but also to spin-parity (1/2)^ states. In other words 
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T]f^ has nonzero overlap with spin 1/2 states. This couphng can be written in 
the most general form as: 

< lip) >= (Ap'^ + B^^Mp) (12) 

Using the condition 7^//^ = one can easily obtain that B = ^^!iZi. There- 
fore, in general case spin 1/2 states give also contribution to the considered 
correlation function but they contribute only to the structures which contain a 
7^ at the beginning or which are proportional to p'^. By choosing the ordering 
7m ^'li^ ^ ^iid ignoring the structures proportional to p'^ and the structures 
that contain a 7^ at the beginning, the contributions to the correlation func- 
tion from states that have spin 1/2 are eliminated. 

Now we turn our attention to the calculation of the correlator function from 
the QCD side for the large Euclidean virtuality p'^ = (p — q)^ << in terms of 
the nucleon distribution amplitudes. Using the explicit expression for the 
isobar interpolating current Eq. (JSj) and axial current Eq. ([6]), and carrying 
out all contractions for the correlation function, Eq. (j4j), in x representation 
we get the following result 

n,.(p, q) = j —^{{C^,U{l.l,)p.e^'' < I [Au^^{M{x)d%{Q) 

-Au''^{Q)v!l{Q)d%{x){2g^rigxega<t>{^)pp + go.vgpega^{^)xp)\ I N{p) >} (13) 
where we have used the light cone expanded light quark propagator as PU] : 



-"^'^-^'ii^] (14) 
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The terms proportional to the gluon strength tensor can give contribution 
to four and five particle distribution functions but they are expected to be 
small PHI25II26] and for this reason we will neglect these amplitudes in further 
analysis. The terms proportional to < qq > can also be omitted because Borel 
transformation eliminates these terms and hence only first term in Eq. f|T4l) is 
relevant for our discussion. From Eq. f[T^ it follows that for the calculation of 
^ij.u{p, q) we need to know the matrix element 

< I [4«^(0)«^(a;)d^(0) | N{p) > . 

This three quark matrix element between vacuum and the proton state is given 
in [211251126] as: 

abc„,a I 



'in 
T 



2 2 
X TTl/ 



2 2 

+ AAm\{:fe'^^C)apN^ + A->m\{'^ ^'^^C)ap{i(T'^'' XyN).i + AQm\{:fe'^^C)ap{i:N).i 

2 2 



where the calligraphic functions are functions of the scalar product {px) and 
of the parameters a^, 2 = 1, 2, 3. These functions can be expressed in terms of 
the nucleon distribution amplitudes with increasing twist. Explicit expressions 
of distribution amplitudes with definite twist are (see also [8ll24f25ll26j ): 

Si = Si, 
2pxS2 = Si- 5*2, 
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P1=P1, 

2pxV2 = Pl-P2 



(16) 



2pxV2 = Vi-V2-V3, 

2V3 = \/3, 

4pxVi = -2Vi + ^3 + ^4 + 2^5, 
4pxV, = V,-Vs, 

iipx^Ve = -V, + V2 + V3 + V^ + V,-Ve (17) 



Ai = Ai, 
2pxA2 = -Ai + A2-A3, 

2^3 = ^3, 

ApxAi = -2Ai - A3 - + 
4pxA5 = ^3 - ^4 

4:{pxfA6 = - + A3 + - As + As (18) 





= Tu 


2pxT2 


= T^ + T2-2T3, 


2T3 


= Tr, 




= T,-T2-2Tr, 


2px% 


= -Ti + T5 + 2T8, 




= 2T2 - 2r3 - 2T4 + 2r5 + 2T7 + 2r8, 


ApxTj 




4.{pxf% 


= -Ti + T2 + T5 - Te + 2T7 + 2T8 



where Eqs. ( fT6l) . (|T71) . (|T8l) and (fT9l) are for scalar, pseudo scalar, vector, 
axial vector and tensor distribution amplitudes respectively. The distribution 
amplitudes F = Si, Pi, Vi, Aj, can be written as: 



where with i = 1, 2, 3 corresponds to the longitudinal momentum fractions 
carried by the quarks. 




(20) 
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Using Eqs. (fT3ll - (l20l) and after carrying out the fourier transformation, the 
correlator function is obtained in the momentum representation. Its exphcit 
expression is given in Appendix B. 

In order to construct sum rules for axial N — A transition form factors we 
need to choose structures. From Eq. f|TT]) it follows that in principle different 
tensor structures can be used for obtaining sum rules for axial form factors. 
We have chosen the structures proportional to g^u ^' hi IfiP'u ^' ■• 9iiv 
QfiQu i>' to obtain sum rules for the form factors C^^ C^, + and C^, 

respectively. 

Choosing the coefficients of the structures g^j^^ f ^, g^p'^, gfj,y f , q^q^ f in 
Eqs. (11 II) and fIB.ll) and applying the Borel transformation with respect to the 
variable p'^ = {p — qY which suppress the contributions of the higher states 
and continuum we get desired sum rules for the form factors C^^ C^, and 
C ^ as! 



C,{Q') = j dx, J dx, "" [2V, - Ti] {x. 



to 



+ J dX3 J dXi'—— Ti(x-)+ j dXs J dXi J -^e ^^^(Xg - ti)T234578(x.) 



X J J J t iv± 

to ^ to to ""^ ^"-^ 

dx-i j dxie ^'b ^ ^2 (~^o + 3^3)^234578 (a^D ) (21) 

to 



2mN 



m 



V3X 



1 1—2:2 
dx2 

to to 

2m7v(l - 2ti) 



^2 



s(ti,Q-') 



dxi I dti e "'b 

A123 + Ti27){Xi) 
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+ 



Am%(-l + ti)(ti-X2) 



1—X2 



+ / dx2 / dxi e 



to 



125678 



)(Vi23-Ti23)(a;,) 



+ 



2mjvto(-l + 2to) 



[-A123 + 7127) {Xi 



{-l+to){-X2+to) 



g2 + m%tl 
tl{4x2 - 5{1 + X2))to + 6tl) 



+ 



(g2 + m%tl) 



1 1-3:3 
+ / dxs 



X3 



s(ti,Q^) 



to 







dxi J dtie ^s" 
to 



(^123456 + ^25678) (a^j) 



■-^(2Vi23 + A23)(a;J ^727 ^i27(a;J 



I mi,(-l + ti)(ti-X3) , 

H TTT^ (,—^1^123456 " ^^123456 + ^ -'125678J l^Tj 



1-13 



+ 

+ 
+ 



j dxs j dxie 



rriNtl 



+ m^^ig 



(2Vi23 + A23)(^D 



to 

/ 2m%{-l + to)4{-xs + to) m%tl{Axs - 5(1 + ^3)^0 + 6tg) 



2mArto(-l + ^0 
Q^ + m%tl '''' 



\ (g2 + m%4y 

mlf{-l + to)to{-X3 + to)' 
MUQ' + m%tl) , 



(Q2 + m%tl)^ 

(2Vi23456 ~ ^123456 + 2Tl25678) (^^i) 



(22) 



n\ ( 1 1-X2 



a/3Aa 



C5(g') = e^B <^ /da;2 / dx^e + _ + A3)(xi) 



to 



1 l-x-3 



+ 



j dxs J dxie 

to 

1 I—X2 X2 



+ / (ia;2 / dxi I dtie ^^^'^ 



to 



to 



ti 



X 



+ 



^ Q^ + m|,(-l + 2ti) + 5(ti,Q^) ^| m^(ti-X2) ^ 

^ — ^^123[Xi) H 770 ^12567812^1; 



Ml 



Ml 



Q^ + m%{-l + 2ti) + s{ti,Q^ 



2M| 



1 1-X2 



+ /cix2 / dx,e-^s. --^0 

to 



g2 + m%tl 



-2^23 + ^123 + ^127) (a^i) 
+ m%to{to - a;2)Ti25678(2;; 
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{-m% + Q^ + s{to, g') + 2m%to){-Vi23 - A23 + ^123 + ^127) (a^i) 

1 1—13 X3 q2j 2 

+ jdx^j dx^j dhe'^ '^''^^' {2V,23M) - A,23Am{x[)) 

to to ^ ^ 



1 1— S3 so 3 /' ^ 



to ^ 

■^=^^|i^c.w^Ma (23) 



'"a 



I-X2 



X2 



«(«i,Q^) 



(ixi / dtie 



to 



to 



94 ^-r V ^ , -2m^(-l + ti)(-l + 2ti) ^ 
2^123 + Ti23)(a;i) H —^-^ ^127(2;^) 



+ 



2m3^(2(-l + ti)% - (1 + 2(-2 + ti)ti)x2) 



M%ti 



Mm 

(^123456 + ^25678) (2^0 



I—X2 



+ 



j dx2 j dxie 

to 

mAr(-l + to) 



4mjv(-l +to)tc 

g2 + ^2^t2 



A23(2;iJ 



+ 
+ 



^ . X 2mjv(-l + to)(-l + 2to) ^ 
g2 + m%tl ^^^^^^^^ g^ + m^tg ^^^^^^^^ 

'8m%{-l + tof4{-X2 + to) 4m^(-l + to)to(3a;2 + to(-4 - 5x2 + 6to)) 



V (g2 + m^t2)3 

4m^(-l + to)2(-x2 + to)\ 



M|(g2 + m^t2) J 

1 1—2^3 a^3 

+ j dxs j dxi j dtie "s 
to to 



(^123456 + ^25678) (^i) 



f^ mjv(l-ti) 



+ 



+ 



2m%{l-h){h-X3) 



1 



M|ti 

1-13 



(Vi23456 — -^123456 + ^25678) {x'i) 



j dxs J dxie 



to 

+ 2m 



^ mjv(-l+to) 



toAl23{x'i)+Ti23ix'i) - (-1 + 2to)ri27(xO 



+ 



2 , 2mif(-l + to)4{-X3 + to) tl{3x3 + to(-4 - 6x3 + 6to)) 

"-AT 

-l + to)(-a;3 + to)' 



(g^ + mM' 



(g^ + m%4 



Ml 



(Vi23456 — 2^123456 + ^25678) (^^i) 



(24) 
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where following [25], the following shorthand notations for various combina- 
tions of the DA's are employed: 

^34578 =T2 — T3 — T4 -\- + Tj + Tg 
^23 = Ti +T2 — 2T3 
7^27 = Ti — T2 — 2T7 
^25678 = + T2 + T5 — Tg + 2T7 + 2Ts 

V,23 = Vi-V2-V3 
Vi23456 = -Vi + 1^2 + ^3 + ^4 + ^5-^6 

^123 = - A2 + A3 
>4i23456 = - A2 + A3 + - A + As (25) 



and 



J^{x.i) = J^{xi,X2, 1 - Xi - X2) 

:F(x-) = J^(xi, 1 - xi - X3,X3) (26) 

The Borel transformation and the continuum subtraction are performed using 
the following substitutions (see plp5] ): 



^ to 

P(t) = f dt Pit) ^ J_ I ^o(t)e~^ + Pi^^)^'^ 

to 

dt P^^) = I ^ P^^) ^ ^ / ^p(t)e"4 
(g - tpf I t^s- p'y 2M% i t3 



p(to)e ^ tl 



2to(g2 + tgm^)M2 2(Q2 + tlm%) 



d p{to 



dto to{Q^ + tlm%] 



s(t,Q') = il-t)M'+^-^-^Q', 



2^ _ V + ^0 - + "^mj^Q^ - (g2 + _ 

2m^ 



.2 



where = — and to is the solution of the s{to, Q"^) = Sq. The terms ~ e 
are so-called surface terms which appear in successive partial integrations 



13 



to reduce the power of denominators. In the hadronic representation for the 
correlator functions the Borel transformation is performed by substituting 




3 Numerical analysis 

From explicit expressions of the sum rules for the axial N — A transition form 
factors it follows that the main input parameters are the nucleon distribution 
amplitudes (DA's). In general these distribution amplitudes contain hadronic 
parameters which should be determined by some means. Various methods to 
determine these parameters give different results. In this work, we considered 
all three different determinations of these parameters: a) QCD sum rules based 
DA's, where corrections to the DA's are taken into account and the parameters 
in DA's are determined from QCD sum rules, b) A model for nucleon DA's 
where parameters are chosen in a such way that the nucleon electromagnetic 
and axial form factors are described well within LCSR and c) Asymptotic 
forms of DA's of all twists (see for example [23] ). 

Explicit expressions of corresponding DA's can be found in [24] and for com- 
pleteness we present their expressions in Appendix. In the appendix, the val- 
ues of the hadronic parameters obtained from three different methods is also 
presented. 

For the numerical evaluation of the sum rules for the N — A transition form 
factors we need also specify the values of the residue of A baryon A a, the 
continuum threshold Sq and Borel parameter Mg. The residue A a and Sq are 
determined from analysis of the mass sum rules: Aa = 0.038 GeV^ and Sq = 
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2.6 — 3 GeV'^ |27f28ll29ll30f31j which we have used in numerical calculations. 

The Borel mass parameter is the auxiliary parameter of sum rules. There- 
fore we need to find a suitable region of M^, where physical results are in- 
dependent of this parameter. A suitable region of is determined in the 
following way. From one side, has to be small enough in order to guaran- 
tee suppression of higher resonances and the continuum contributions to the 
correlation function and from other side it should be large enough in order 
to guarantee convergence of the light cone expansion with increasing twist in 
QCD calculation. 

In Figs. [H [31 andlU we present the dependence of the form factors C^^Q"^), 
CiiQ^), C^iQ'^) and Ce{Q'^) on M| for first set of DA's with two fixed values 
of So and three fixed values of respectively. From these figures it follows 
that all considered form factors exhibit good stability with respect to variation 
of in the region 1.2 GeV'^ < < 2 GeV'^, so this region of can be 
considered as a working region where form factors are practically independent 
of M^. We performed similar analysis for the two other sets of DA's and 
obtained that within the above mentioned region of M^, the form factors are 
rather stable to variation of M^. 

In Figs. El El El and [S] we present the dependence of the form factors 
i = 3,4,5,6 on at fixed values of and sq for all three sets of DA's. 
From these figures we see that form factors are very sensitive to the choice 
of DA's. For form factor C^^IQ"^) sets 2 and 3 leads to the same result, for 
Ci{Q'^) and C^IQ'^) form factors sets 1 and 2 give results which are very close 
to each other and for CqIQ"^) up to = 4 GeV^ all sets of DA's leads to 
the different results and when > A GeV^ all three sets of DA's leads to 
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indistinguishable predictions. Our results on form factors C4(Q^) and C^^Q"^) 
for three sets of DA's satisfy the relation C^iQ'^) = C^iQ'^)/^ assumed in the 
experimental analysis. 

Figs. ([7]) and ^ contain also the predictions of lattice calculations from [T3]. 
The points correspond to the values obtained using a hybrid action and assum- 
ing Za = 1-1 (see pj] for details of the notation). Note that, due to different 



conventions used, the lattice results has been multiplied with y2/3. From the 
figures, it is seen that there is a good agreement between the lattice results 
and our predictions within error bars. Note also that, the biggest difference 
between sum rules predictions and lattice calculations is seen in C^i^Q^) when 
the asymptotic distributions are used, i.e. when the hadronic parameters from 
Set 3 are used. 

And finally, in Fig. ([9]), i?(Q^)is plotted as a function of Q^. The function 
R{Q'^) is defined as: 



From Eq. ([3]), it is seen that, assuming PCAC and pion dominance, R{Q^) = 1. 
From Fig. ([9]), it is seen that at ~ 1.5 GeV, PCAC and pion dominance 
approximations are valid. But for larger values of Q^, R{Q'^) deviates consider- 
able from unity, and hence we may conclude that PCAC and pion dominance 
assumptions break down. 

As we have already noted, the axial form factors for the A ^ N transition 
have also been calculated in the framework of the constituent quark model in 
[Ti] . Our predictions on these form factors differ from the results of [14j . 



(28) 
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In summary, we have calculated axial — A transition form factors which 
play important role for understanding the dynamics of weak A'^ — A transition 
and compare our results with existing lattice and quark model predictions. 
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A Appendix A 



In Eqs. ffTB]) . ffTTI) . ffTSl) and ffT^ the distribution amplitudes depend on scale 
and can be expanded with the conformal operators, to the next-to-leading 
conformal spin accuracy, they are obtained in p^25ll26] : 







/i) 


= 12OXiX2X3[0°(/i) + - 80:3)], 


V2 




/i) 


= 24xiX2[0^(/i) + 0^(/i)(l - 5x3)], 






n\ 
fl) 


— -L^a^31^4l/^Jl-'- + r4 l/^JFl + — X3(^i — X3JJ 








^'xb't — x-x — 10x1X9)1, 




'x- 


u) 


= 3i^?fu)fl — x-^) + '0rfu)[2a:iX9 — x-^il — x-^)] 








+V^5+(/i)[l-X3-2(x?+x2)]}, 






/i) 


= 6X3[0°(/i) + 05+(/i)(l-2x3)], 






/i) 


= 2[0O(/x) + 0+(/i)(l-3x3)], 


Ai 


^Xi 


/i) 


= 12OXiX2X303 (;U)(X2 - Xi), 


^2 




/ij 


— Z4l:XiX2<P4 l^/iJl^X2 — XiJ, 


^3 




/i) 


= 12X3(X2 - + ^+(/i)) + ^4~(/^)(l - 20:3)}, 


^4 




/i) 


= 3(X2 - Xi){-V^5(/i) + V^^(/i)x3 + V^^(/i)(l - 2X3)}, 


A, 


^Xi 


/i) 


= 6X3(X2 - Xi)05 (yU) 


A, 


yXi 


/i) 


= 2(X2 - Xi)0g (/i). 


Ti 


^Xi 


/i) 


= 12OXiX2X3[0!i(/i) + -(03 - 03^)(/")(l - 3X3)], 


T2 




/i) 


= 24xiX2[e4°(/i)+e4''(/")(l-5a;3)], 


T3 




/i) 


= 6x3{(e4 + 04 + V^4)(Ai)(l - 2:3) + (C + 0r - V^4")(/^)[a;? + a;2 - 2:3(1 - 0:3)] 








+ {^t + 01 + V^4^)(/U)(l - 2:3 - 10X1X2)}, 


T4 


^Xi 


/i) 


= ^{(^5° + 05 + V^5)(/^)(l - ^3) + (4" + 05 - V^5")(/^)[2a:iX2 - X3(l - X3)] 








+ (a + ^t+ ^5^)(/")(l - ^3 - 2(X? + X^))}, 


T5 




/i) 


= 6X3[e5°(/i)+e5+(/i)(l-2X3)], 


Te 


yXi 


/i) 


= 2[0°(/i) + i(06 -0+)(/i)(l-3X3)], 


T7 


^Xi 


/i) 


= 6X3{(-e4 + 04 + ^4) (/")(! - 2^3) + (-C + 0r - ^4")(/")k? + ^2 - 2^3(1 - 2:3)] 








+ (-^4^ + 01 + ^^)(/W)(l - X3 - 10X1X2)}, 


Ts 






= ^{(-^5° + 05 + ^5)(/^)(l - ^3) + (-4" + 05- - i^^){^l)[2x^X, - X3(l - X3)] 



^1 



+ (-^5+ + 05^ + ^5^)(/")(l - 2:3 - 2(X? + X^))}, 
/i) = 6X3(X2 - Xi) [(^4° + 0^ + i^l + ^4^ + 01 + V^4^)(/i) + (C + 01 " ^4-)(/")(l " 2^3) 
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S^{X,, /X) = %2- Xi) [- + <Pl + (/i) + + 05 - ^5 



2^3 



05^ - < 



2x3) 



(A.l) 



the parameters used above are defined in terms of the following eight inde- 
pendent parameters /at, Ai, A2, V^, A", /J, /J and as 



03 = 06 = /a^ 

04 = 05 = ^(Ai + /iv) 

^0 — ^0 - "'"X 

?4 — ?5 — g^2 
^4° = ^5°=^(/iV-Al) 

03 = f ^1 , 

03^ = ^(1-3^'^), 

0r = ^ (ai(i - vt - m + fN{'2A\ - 1)) , 

01 = ^ (Ai(3 - 10/f) - fr,{Wt - 3)) , 

= (Ai(2 - 7ft + /r) + /iv(A« + Wt - 2)) , 
i^t = (Ai(-2 + bft + 5/r) + /iv(2 + bA\ - Wt)) , 
C = ^A2(4-15/2^), 
e4^ = ^A2(4-15/2^), 
05=^(Ai(A'-/r) + /iv(2A«-l)) , 
05+ = (Ai(4/f - 1) + /^(3 + 4^,'^)) , 
V'5~ = ^ (Ai(/i' - /r) + /iv(2 - - ^Vf)) , 
i^t = -| (Ai(-1 + 2ft + 2/r) + /Ar(5 + 2^5* - 2Vt)) , 
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— "^-^2/2 ) 

</'6 = I (Ai(1 - 4/f - 2/r) + /^(l + 4A^)) , 

<l>t = - Vt) + /iv(4n' - 1)) (A.2) 

Our numerical values are obtained using: 

fN = (5.0 ± 0.5) X 10-^ GeV^, 
Ai = -(2.7 ± 0.9) X 10-^ GeV^, 

A2 = (5.4 ± 1.9) X 10-^ GeV^ (A.3) 



And for other five independent parameter we have used three sets as: 



Set 1 

a;' = 0.38 ±0.15, 
V;'' = 0.23 ±0.03, 
/f = 0.40 ±0.05, 
/2^ = 0.22 ±0.05, 

= 0.07 ± 0.05 (A.4) 



Set 2 



1 

"14' 
13 

~42' 
/f = 0.40 ±0.05, 

/2^ = 0.22 ±0.05, 

/i" = 0.07 ± 0.05 (A.5) 



Set 3 



0, 
1 



3' 
4 



15' 
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/r = ^ (A.6) 

Note that the asymptotic forms of DA's can be obtained from Eqs. (1A.1I) using 
the values given in the 3'''^ set. 



B Appendix B 

In this appendix, we present the exphcit form of the correlation function con- 
taining all the Dirac structures. 



1— xi 



mN [/gt,uX2 + iluQiM- ^fl'MZ^)(l - 3:2)] {Si - P1 + V3- As){Xi 



+ 
+ 



2gM(P^ + 2'?i.)(-l + a;2)+fi'Mi^(p'-g + P-('?-a;2j5) + g - (p' + g) )x2^1t>v%^i ^i( 



+ Iv ^Qi^i-l + X2)- i>lv%{^ + X2) - gf,^X2{2p'.q + q^) 
"^qM + 2gi.)(-l + a;2)a;2 + gf,u{p' + g)^a;2} Ti{xi) 
1 



1 l-xi 

+ J dxi j dxs 



{q - xspf 



A3 



mN \i>g^vX3 + {ivq^- igp.v)i^ - x^)] [V^ - —j (x'i) 
+ g/. [(7. ^- 2g,)(-l +X3)- ^7.^3] (^1 + y) ^^^) 

+ [i/lu + 2q^)qf, - gf,^{^' 4+q^)- 2q^{p'^ + g^)x3 + g^u{p' + qfx^ Ti(x-^ 

1 1— Zl X2 

+ j dxi j dx2 j dti 



5m 



{q - tipf 

-1 + - X2)g;, [tip'^ + (-1 + ti)g^] [h / + (-1 + h) ^ Ai2Zi^Q{xi) 



+ Am\{-l + ti)(ti -X2)q^ 

X {p.{q - tip) [ti /7, - (-1 + ti)7, ^ 

+tip', [2 f ^ + p.{q- hp) + + (-3 + 2ti)p'.q + (-1 + t^g^" 

-(1 - ti)q, [2 / ^+3p.{q - hp) + tip'^ + (-3 + 2ti)p'.q + (-1 + ti)q^] } T234578{x'i 
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- 8m^(-l + h){ti - X2)q^, [hp', - (1 - h)qu] [h / - (1 - ^i) ^ T^2^&7M) 



2m^(l - ti){tx - X2)q^i^vAi2Mm{xi) 



(q - hpf 

+ m^(ti - X2) [g^,, (p.{q - hp) + hp'^ + (-1 + 2h)p'.q + (-1 + h)q^) 



234578 



+ {(3 - 2h) /7. + 2(-l + ii) (7. ^ - 2g.) - 2p'J' 

- m%{h - X2) [g,j,v / + (-1 + h){g^,w 4 + Qnlv)] ^125678(2^0 

+ m^g^(-l + h) [1. ^(-1 + h)- illvh + 2p',ii] - + Pi - P2)(a:i) 
+ mjv [2 /7, iq^ - 2 i{gi,^p.{q - hp) -pM(-l + h) 

-2 ( /gf,uP-{q - hp)+ i>lv iq^i^ T^'p^Qt}) h 

+7,g^(-l + h){-^p'.q + 3p.(g - tip) -q^ + {p' + 9)=^^!)] Vi23(a;i) 

- lrn% - ii) [37, ^{1 - h)+ Ml -h) + 2 il^M 

+ [-2g^(l - h) (g.(l - ti) + 2pUi) + - 2(p'.g + + + } (V^4 - Vz){x,) 

- [ {2 + ii) + 2 /ii} {g^,,p.{q - hp) - q^ + 2?.) + 2g^(p', + 

- 7i^gM(-l + ^i) {-P'-? - ^iP) - + (P' + ^23(2^0 

- 2g^M2(-l + ii) + ti) A345(a:;i) 

- g^M2(-l + h) i-2q,) (-1 + h)- ^lvh\ {A-i - A^){x,) 



+ mN 



{4gtiu - luq^?)P■{q - hp){-l + ^i) 



+ / {-2qt,qv + g^vP-iq - hp)h + 2q^{p'^ + q,)h)\ T^i2z{xi) 



+ M 



/Ti^ - ^i) + 2 ilqix^v^ / [gp.vP-{q - hp) - 2q^{p'u + 3?,)) 



+4 ^'qM + 5-)^?+ ^(-1 + ^1) {9^.uP■{q - hp) - 2qM + 2?.) + + q^)ti) 

- luq^^{-^ + h) {p-{q - hp) - + 2p'.q + {p' - q)\) 1 Ti27{xi) 



+ g^m^(-l + h) (7. g^- 4g,)(-l + ii)- /7,ii - 2^^^! ri58(xi) 



+ 



(-(/+ ^^M- - 7.?m)(-2A23 + T123 + Ti27)(Xi) 



2iq - hp)^ 

+ 3miv^^z.(l^4 - (xj) - 2((/+ gOs'^i. - li^qiJ.)Vi23{xi) 



1 1— xi 2:3 

+ J dxi J dxs J dh 



{q - hp)"" 
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2q^ml{l - t,) [q, - {p', + {2 [/ ^ - p'-q + {p' + q)-{q - hp)] %u57s{x'i) 

+ rUNy- (/+ ^^l] (2Vi23456 - ^23456 + 2Ti25678) {Xj)] {h - X3) 



+ 



{q - tipy 



+ m]^{ti - X3) \g^^ (^i> i- p'-q + (p' + q)-{q - hp)) 

- q^inu 4-2q^){l - ti) - q^^-f^ /ii] T2Z4:-,7s{.x'i) - m%{ti - X3){-1 + tl)g^7^Ti25678 

+ 2mNq^,Ui-l + h)+ +ti) +p',ti)Vi23(a;0 

+ lqi^mli-1 + h) [(^7, + 2q,){-l + h)- /7,ti] Vi345{x'i) 

- Iq,ml{-1 + h) [(g^7. + 2g.)(-l + ti) - (/7. - ¥.)ti] {V, - Vs){x'^ 

-5'M<^"^ir(#(-l + -a:3)Vi23456(a^i) 

+ "^jvg^ [^(-1 + h)+ /ti] + ti) + p'^ti] A23(a;-) 

- \q^^rn%{-l + h) ((7, ^- 2g,)(-l + ii)- /7,ii) ^msl^;-) 

- \q^rn%i-l + h) ((7, ^+ 2g,)(-l + ii) - (/7, - Ap'Jh) (A^ - A4)(xD 

- ^mNq^.{-l + h) (/ g^- 2 /q, + 7, [p'.g - (p' + - tip)]) T^six'^) 

+ ^mNq,,{-l + h) i>lv i-2{i) + 2 (i)q^ lv\~ V -q + (p + q)-{q - hp) 



+ 4(/+ + ri27(a;:) - 2q^ml{-l + + h) + p'ji)r,,s{x'^ 



where p' — p — q- The short hand notations for the functions Ax-, Tx and Vx 
are defined within the main body of the text. The shorthand notations that 
do not appear within the main text are defined as: 

^1345 = 2Ai + A3 + - 2^5 

Tiss = -Ti + T5 + 2T8 

V1345 = -2Vx + V^3 + V^4 + 2^5 (B.2) 
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Figure Captions 

Fig. 1 The dependence of the form factor CslQ'^) on the Borel 
parameter squared for the values of the continuum 
threshold sq = 2.6 GeV^ and sq = 3.0 GeV^ and at the values 
of g2 = 3 ± 1 GeV^ 

Fig. 2 The same as Fig. ([T]) but for the form factor C4(Q^). 

Fig. 3 The same as Fig. ([T]) but for the form factor C5(Q^). 

Fig. 4 The same as Fig. ([T]) but for the form factor CqIQ"^). 

Fig. 5 The dependence of the form factor C^IQ"^) for three different 
sets of distribution amplitudes at the continuum threshold 
So = 2.6 GeV"^ and the Borel parameter M| = 1.5 GeV"^ 

Fig. 6 The same as Fig. (jSj) but for the form factor G4{Q'^) 

Fig. 7 The same as Fig. ([5]) but for the form factor G^IQ'^). Results 
from the lattice are also shown. 

Fig. 8 The same as Fig. ([7j) but for the form factor Gq{Q^) 

Fig. 9 The dependence of R{Q'^) as a function of for the three 
sets of DA's. 
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Fig. 1. The dependence of the form factor C3((5^) on the Borel parameter squared 
for the values of the continuum threshold sq = 2.6 GeV'^ and sq = 3.0 GeV^ 
and at the values of = 3 ± 1 GeV"^ 
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Fig. 2. The same as Fig. ([T]) but for the form factor Ci{Q'^). 
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3. The same as Fig. ([T]) but for the form factor C^{Q'^ 
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4. The same as Fig. ([T|) but for the form factor Cq{Q'^ 
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Fig. 5. The dependence of the form factor C3(Q^) for three different sets of distribu- 
tion amphtudes at the continuum threshold sq = 2.6 GeV'^ and the Borel parameter 



M2 = 1.5 GeV^ 




29 



0,6 

0,5 

0,4 

'So,3 
U 



0,2- 
0,1- 
Q 







<> C5 ^ 



— Setl 
■■ ■ Set2 
-- Set 3 



O m^ = 0.594(1) GeV 



0.498(3) GeV 



O = 0.357(2) GeV 




Fig. 7. The same as Fig. ([5]) but for the form factor C5((5^). Results from the lattice 
are also shown. 
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Fig. 8. The same as Fig. ([7]) but for the form factor Cq{Q'^ 
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